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We study the magnetic and weak magnetic dipole moments of charged leptons in the framework 
of the Two Higgs Doublet Model type III (THDM-III) with four zero textures. We first analyze 
the possibility that the discrepancy between the experimental measurements and the theoretical 
prediction of the muon anomalous magnetic dipole moment (a M ) in the Standard Model (SM) can 
be explained in the context of the model and then find values for the parameters of the THDM-III: 
neutral and charged Higgs boson masses {mu, m H ±), mixing angles (a,/3), CP-violating angles 
0i, 62 and an additional parameter arising from the mass matrix (7/), for which a correspondence 
between theory and experiment is possible. We use a'^ HDAI ~ 111 and current low energy processes 
K — K mixing, r —> /qqx, f-i —> eee, t — > /.17, t —> ey, /r —> e-y, b —¥ s-y , B — > D(D*)tu 

and the rare higgs decay h —> t[i in order to determine the allowed parameter space of THDM-III 
and from these we give a prediction for the magnetic and weak magnetic dipole moment of charged 
leptons, in particular of the t lepton. The obtained magnetic dipole moment of the t lepton is of 
the order of O(10 _s — 10 -7 ) and the weak magnetic is of the order 0(1O -10 — 1CF 7 ). 

PACS numbers: Valid PACS appear here 


I. INTRODUCTION 

The study of dipole moments in particle physics has become an important tool both theoretically and experimentally 
in the search of new physics. Motivated in part by the fact that within the SM the electric dipole moment at the one 
loop level is zero, albeit different Beyond the Standard Model (BSM) theories predict a small yet different than zero 
value, this would serve as a new source of CP violation. In regard to the magnetic dipole moments, the theoretical 
prediction of the lias been found to differ from the experimental value by [I. 

Aa^ = a^ xp - a^ M = 288(63)(49) x HT 11 , (1) 

which is greater than 3 standard deviations. This discrepancy between theory and experiment could give a hint to 
new physics. Of the several BSM models that have been used to try and offer a possible explanation to A a^, one 
of the simplest is the THDM, as was done in the work of 0 and [ 3 ], which focused on THDM’s of the type I and 
type II. What differentiates each version of THDM is the way in which the Higgs doublets couple to fermions; in the 
type I only one of the doublets couples to all fermions while the other remains inert, type II has one doublet coupling 
exclusively to the up type quarks while the second doublet couples to down type quarks and leptons. Finally, in the 
type III (THDM-III) both doublets couple to all fermions, up and down quarks and leptons. As far as the authors 
know, no one has determined if the THDM-III is a viable model to explain Aand in turn how it might constrain 
the parameter space of the model. 

Unlike the fi, the r lepton has a short life-time of (290.3 ± 0.5) x 10 -15 s that does not allow measurement of its 
anomalous moments by the spin precession method used with the electron and muon. The dipole moment of the t, 
with a 95% C.L., is [ 3 ] 


- 0.052 < a T < 0.013, 

(2) 

with a central value 


a T = -0.018(17), 

(3) 

and the SM theoretical value is given by 0 


a SM = 1177 _21(5) x 10" 6 , 

(4) 


different BSM models give predictions for the magnetic moment of the r between the range of 10 9 to 10 - 6 0 , 173 , 0 . 
On the other hand, with respect to the weak properties of the r, the prediction of the SM was calculated in [5\ 
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whose value was found to be: 


a ^- SM = _( 2 . 10 + 0i61 j) xl0 -6 i (5) 

different BSM models predict values between the range of 10 -10 to 10 - 6 mmm- The real and imaginary part of 
the Weak Magnetic Dipole Moment (WMDM) of the r, with a 95% C.L., were reported in [IT] and are given by 

Re (cl™) < 1.14 x 10“ 3 , Im (a™) < 2.65 x 1(T 3 , (6) 

while the Weak Electric Dipole Moment (WEDM) is 

Re (d]T) < 0.5 x 10 _17 ecm, Im ) < 1.1 x 10 ~ 17 ecm. (7) 


The ffV effective vertex describing the interaction between a neutral vector boson and two on-shell fermions can be 
written as: 

r f J v (q 2 ) = ie [ 7 „ [F$ - Fy 75 ) - + E^ 75 ) + oprf ( iF % + Fl 7s )] , (8) 

where q = p 2 — p\ is the four-momentum. The dipole moments are defined as follows: 


and the weak dipole moments 


a f = d f = -ei^(O), 


7 = - 2 m f FZr z (m' z ),dy = -ei^= Z (”4)- 


(9) 


( 10 ) 


In this work, we focus on the contributions to magnetic and weak magnetic dipole moments of charged leptons that 
arise from a THDM-III with a four zero texture. In particular we study the anomalous magnetic dipole moment of 
the t which can be measured through the decay r~ —> l~ v T VF 1 m and its weak magnetic dipole moment m 
The structure of the paper is as follows; in section [TT] we introduce the Lagrangian for the Yukawa sector of the 
THDM-III and how it can be expressed within the four zero texture framework. Afterwards in section [YTT] we present 
our calculations for the electromagnetic and weak dipole moments along with there one loop Feynman diagrams. While 


in section IV we present the effective Hamiltonians and the current values on experimental bounds and measurements 


for the Low energy processes considered. Results and numerical analysis are given in section[V] Finally, our conclusions 
are presented in section |VT 


II. THE TWO HIGGS DOUBLET MODEL WITH A FOUR ZERO TEXTURE 

In this section we obtain the coupling of the Higgs bosons to fermions H(H ± )f i fj. We show the way in which we 
diagonalize the fermion mass matrix, where the structure of the Yukawa matrices is such that after summing them 
we recover a mass matrix with a 4 zero texture. 

The general THDM, has been studied under different theoretical frameworks, such as the Minimal Flavor violating 
(MFV) THDM P3, which provides the minimal level of FCNC consistent with data. The THDM with Alignment is 
a scenario in which flavor violation does not occur at tree level; the contributions coming from vertices where flavor 
change occurs are not present and therefore could be smaller then the ones present in our model. 

The assumption of textures implies a specific pattern of FCNC Higgs-fermion couplings, known as the Cheng-Sher 
ansatz 'I rnirn :j . These vertices satisfy FCNC bounds for Higgs masses lighter than O(TeV), which is one of the reasons 
why we work in this formalism. Secondly, it gives the structure of the Yukawa matrices as a function of fermion masses 
in such a way that the elements of the CKM matrix are reproduced. But there are different possibilities for the texture 
forms, such as a textured six, five or four zero mass matrix which are ideas that has been extensively developed; one 
of the first efforts to work in this scheme assumed a six zero texture type [T5], nonetheless this texture has been 
ruled out 1 1 16]. In this work we focus on the four zero type motivated by the fact that they currently satisfy many 
phenomenological observations. Furthermore, for certain four zero texture types as the Nearest Neighbor Interaction 
(NNI) is obtained through a Z 4 flavor symmetry in [171 . For a more in depth study of the formalism of textures and 
the THDM, we suggest Refs. [TS], [IS], [2D] . 

The Yukawa Lagrangian of the THDM-III given by: 
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4 = Y?Q l * iUr + Y 2 u Q l $ 2 u r (11) 

+ Y^Q L ^id R + L^^dR 

+ YiLr^iIr + Y^Li,^21r + h.c., 

here $; = (4, <74 denote the Higgs doublets, Y/ (/ = u,d,l) are 3x3 Yukawa matrices, = ia 2 ^j and 
Q l = (ul, dz,). Since both Higgs doublets couple to all fermions, the mass matrix for each fermion type / receives 
contributions from both vevs v± and V 2 , i.e. 


Mf = 71 ( VlY * + V2Y 2 f ) ■ 


( 12 ) 


In order to obtain physical fermion masses, we need to diagonalize the mass matrix; this is achieved through a 
bi-unitary transformation which is given by: 


M d = UfMfU) = U f ^~ (v!Y( + v 2 y[) 


72 

where the diagonalization matrix that we are using is obtained from El 


ul 


(13) 


U f = 0}P f , 


(14) 


in which Of is of the form 


°/ = 


(J 
V 

\s/A 


m 2 m 3 (A — m\) 


A(rri2—mi)(m 3 — mi) 


mi (mi —A) 
(m 2 —m 1 )(m 3 -m 1 ) 
(A-m 2 )(A-m 3 ) 
A(m 2 —mi)(m 3 —mi) 


/ 

-/ 


mim 3 (m 2 - A) 
A(m2-mi)(m 3 -m 2 ) 


m i m 2 (A —m 3 ) 
A(m 3 — mi)(m 3 — m 2 ) 


\ 


m 2 (A —m 2 ) 
(m 2 -mi)(m 3 —m 2 ) 
m 2 (A—mi )(m 3 — A) 
A(m 2 —mi)(m 3 —m 2 ) 


/_ 

/ m 3 (A-mi)(A-m 2 ) 

V -^( m 3 m l)( m 3 m 2) / 


m 3 (m 3 — A) 

(m 3 —mi)(m 3 —m 2 ) 


(15) 


such that A. = 70 . 3(1 — 7 /^ 2 ) with f = l, u, d and r 2 = ((p, and 


1 0 0 
P = | 0 e Wl 0 
0 0 e 1 ® 2 


(16) 


The neutral Yukawa Lagrangian given in Eq. (11) can be expressed in terms of the mass eigenstates h°, H°, A 0 as 


tr = 4—V 

2 \m w J 

+ 4~ V 

2 \ m w J 

+ 

2 \m w J 

+ 4—1® 

2 \m w J 


cos a r 

-c 

cos /3 


6 U , + ("AY) (yl) 

g cos (3 \ mi / V 2 / IV 


sin a 
sin (3 


\ mi J 


l'H° 


COS P 

-tan (35 w 


\/2cos(q - 0) / mw\ ApA 
gcos/3 \ mi ) \ 2 ) IV 


V m l J 


^ f (yl 

g COS/3 v TTH J 


’ W 

7 5 l'A° 


I'hP 


( 17 ) 


S m , + ^ sin(a ~fl ( ^ (yA 

g sin p \ m u J \ / uu' 


u'H° 
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9 (rn^\_ Lcoea + V2co S ( a -(3) /m*A N 

2 \mw) sin/3 gsin/3 \ m u ) ' -'nit' 

* f J- cat/« w + (Ay) (yA 

2 ym^/ / gsin/3 \ m ^ / V A )«»' 


«'/i u 


7 VT 0 , 


where F/ = ^ M p — ^F/ and F/ = UfY/Uj. The type-down quarks part is similar to the lepton part with the 
exchange l —> d and mi —»• m,d- 

On the other hand, the charged Yukawa Lagrangian is given by: 


^ch 

•'Y 


= 

di (y^ sin(3 + Y 2 cos/: 

i) U jH 

+ ui (Y 2 cos/3 — Yisin/3 

) djH+ 

+ 

di (y^cosP — Y^sinji 

i) Uj H- 

+ ui (Y^sinf. 3 — Y 2 cos/3 

) djH+ 


Pr 


Pl- 


(18) 


We will use a notation in which r^j corresponds to the couplings ffH in the Lagrangians, for example: 


*r -f (£) 


COS 


cos/3 


a + VijMa-P) /TOjA , ■> 

0 gcosp \ mi 1 \ nv 


\ m 


(19) 


In this work, we assume that the mass matrix has a four zero texture form 


so 


where 


/ 0 D f 0 \ 
M f = D) C f B f , 

V 0 B f A f J 


( 0 Df 0 
D} C f B f 
\ 0 B* f A 


Vl 

V2 


0 dx 0 
d\ ci b± 
0 bl cii 



0 d 2 0 \ 

d* 2 c 2 b 2 , 

0 62 °2 / 


/ 0 di 0 \ 

Ci = d* Ci bi I , i = 1, 2. 
V 0 b* a t ) 


To find the expression for the mass eigenstates, we use the following matrix invariants: 


( 20 ) 


( 21 ) 


Tr(Mf) 

= Cf + A f = m h + m h + m f% , 

(22) 

det(Mf) 

= -D 2 f A f =m h m f2 m h , 

(23) 

x(M f ) 

= B f + CfA f — Dj = m fl m h + m fl m h + m/ 2 m/ 3 , 

(24) 


from these expressions we find a relation between the components of the 4-texture mass matrix and the physical 
fermion masses 
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Af 

Bf 

Cf 

Df 


m h (1 - r-27/), 


m h\ 


I r 2 'yf (r 2 7/ + n - 1) (r 2 7/ + r 2 - 1) 


1 - r 2 7/ 
TO /s ( r 27/ + J*i + r 2 ), 
m h m h 


1 - 7-27/ ’ 


(25) 

(26) 

(27) 

(28) 


where m,f i is the fermion mass and the 77 stands for 77 = 

So far, we have found a relation between the textured mass matrix and the diagonal physical mass matrix, but the 
coupling of fermions to the scalar is given through the Yukawas, whose elements Oi >2 , 6 p 2 , Ci j2 and dy 2 can not be 
simultaneously determined. But instead of having them as free parameters we can do the following: 

Case 1) We will assume that the elements of Y\ are negligible and so 


/ 0 D f 0 
M f = D* f C f B f 
V 0 B) A 


V2 

V2 


0 d 2 0 \ 

d* c 2 b 2 

0 b 2 a 2 J 


(29) 


Case 2a) Following [30], we will consider the case in which each Yukawa contributes to only one element of the mass 
matrix 


/ 0 D f 0 
Mf = ( D*f Cf B f 
V 0 B*f A 


Vl 

V2 


0 d 1 0 

d\ ci bi 

0 b\ 0 


V2_ 

V2 


0 0 0 \ 

0 0 0 , 
0 0 a 2 J 


(30) 


or 

Case 2b) 


/ 0 D f 0 
Df C f B f 
V 0 B} A 


Vl 

V2 


0 0 0 
0 Ci 0 
0 0 0 


V2 

V2 


0 d 2 0 \ 
d* 2 0 b 2 \ 
0 b 2 a 2 J 


(31) 


(u) 


//' 


in both cases, the structure of the Yukawa matrices textures is such that after summing both we recover a mass 
matrix with a 4 zero texture. In working under any of these assumptions we avoid introducing new free parameters 

instead of explicitly determining the Yukawa matrix elements as functions of the fermion masses and 7 /. The 

for the paralell case is shown in the App VIIA and for cases 2a and 2b are given in [30]. 

In regard to the Renormalization Group Equations (RGE) stability, a study of the RGE in the general THDM-III 
has already been performed by the authors of EH, they found a remarkable stability of the suppressed FCNC Yukawa 
coupling parameters. Furthermore, in [32] they explored if four zero textures are stable under the running of energy 
scales and found that the texture zeros of mass matrices Mf are stable against the evolution of energy scales. 


III. DIPOLE MOMENTS IN THE THDM-III WITH FOUR ZERO TEXTURE 

In this section we give our expressions for the magnetic and weak magnetic dipole moments of charged leptons 
within the THDM-III. 
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A. MUON ANOMALOUS MAGNETIC DIPOLE MOMENT 

The Feynman diagrams that contribute to the magnetic dipole moment of charged leptons in the framework of the 
THDM-III are shown in FIG. [I], 



Figure 1: Feynman diagrams that contribute to the magnetic dipole moment of charged leptons in the framework of the 
THDM-III. In the case of the contributions coming of Barr-Zee diagrams we only consider the dominant contribution 
given by diagram (c). 


where li = e, y, r, and lj refers to a lepton of different flavor, H = h°,H°,A° are the neutral light, heavy and 
pseudo scalar Higgs bosons, respectively. The vertices coupling Higgs bosons to charged fermions are obtained from 
the Lagrangian in eqs.( l7p8 l, except for the vertices coupling Higgs bosons to gauge bosons and Z to neutrinos [55] . 
The contributions to the magnetic dipole moments of leptons are CP-conserving, therefore the phases 6\ 2 are set to 
zero. 

We perform the calculation in the unitary gauge via the Feynman parameters. The one-loop contribution is given 
by: 


n PlP2P3 _ 
a li 


E 


1= 1,2,3 

1 '=2 


^PlP2p3 

'll' 


2 

m k 


V&TT 2 


1 1-x 

j dx J dyG p k lP2P3 {x , y ), 
0 0 


(32) 


where the label p n refers to the particles circulating in each diagram, and the index IV denotes the entry of the 
Yukawa interaction, in the case for the aV = 2, i.e. 12 with l = 1, 2, 3. Notice that our results can be ap¬ 
plied to any of the charged lepton. The label V P ) P2P3 is the coupling obtained from the neutral Lagrangian term of UH. 

The G k lP2p3 (x, y) function for the diagram (a) is given by: 


Ga ljlj (x, y) = {x + y)(m;. - m h (x + y- 1))/M %, 


(33) 


with M% = —m%(x+y—l) + (x+y)(m‘f. +mf.(x + y—l). Here lj = e, y, r. For the Feynman diagram (b): V P1 , P2P3 = 1 
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and G^ lP2p3 (x) function is given by: 


Gb ±H±v {x) = 2 m^x/M 2 , 


(34) 


with M 2 = ( mf.x — m 2 H ±). The THDM also has Barr-Zee two-loop contributions to the muon anomalous magnetic 
moment. The dominant contribution is given by [24j 


two loop _ 

a /l 


o ? 4° 

a 2 

8t r 2 s^ rriyy 


N lQ') r ff( r f) r ifp 

f-t, T, b 


(35) 


where rj = (toj/to^), rrif is the fermion mass, N£ = 1(3) for leptons (quarks), Qf is the electric charge of fermions 
and is given by the Lagrangian coupling the CP-odd scalar to fermions: 


C. = h 5 f A 0 , 

2mw 


and finally 


/O) 



l °9(^Y=y)) 

x-y(l- x) 


dy. 


(36) 


(37) 


Because of the term r/ = (mf/m^) there is an inverse mass dependence with the pseudo-scalar, which in our 
work was evaluated in the range of 200 — 1000 GeV, hence the contribution of this diagram is of the order (P(10~ 13 ). 
In a previous study by .25], the Barr-Zee contributions are considerable because they consider pseudo-scalar masses 
between 10 -2 — 10 2 GeV. However, a collider search for the neutral Higgs bosons in the framework of the THDM’s 
rules out little masses for to^o. An OPAL analysis [25] using the LEP-II data up to yjs = 189 GeV exclude the region 
12 < rriA < 56 GeV. 

Therefore, the total contribution to in the THDM-III is 


THDM—III _ H ± H ± v , n Hl J l i _l n twoloop 

~ a p ~r / y 

H—h °, H °, A 0 

lj=e, fi, t 


(38) 


B. Magnetic and weak magnetic dipole moments of the r 

The analysis that follows can also be applied to any of the charged leptons k = e, y, r. Nonetheless, we will 
concentrate on the case Zj = r. The Feynman diagrams that contribute to the magnetic dipole moment of the r are 
similar to the ones used for the muon, except that we now have that m/ i = m T and V = 3 because the entries of 
V P ] P2P3 are 13, i.e. 13, 23, 33. 

For the weak magnetic dipole moment, we also use diagrams like those shown in FIG. [l] with the photon replaced 
by the Z boson, but we must also consider the diagrams shown in FIG. [2] 

The contributions to WMDM of the charged leptons of each Feynman diagram are given by: 


1 1 — X 

( a u) PlP2P3 = X! A a,b,c,e J dx J dyG p k lP2P3 (x,y), (39) 

a, b, c, e q q 

again the label p n refers to the particles circulating in a particular one-loop diagram, A ai Cj e represent terms 
associated to each Feynman diagram, given by 









Figure 2: Additionally to the Feynman diagrams (a) y (b) there are contributions from (d), (e) y (f) to the weak 
dipole moment of the charged leptons a) v . In the Feynman diagrams (e) H = h°, H°, A° and (f) H = h°, H°. 


Ab 


_.PlP2P3 

'III' 


47 T 2 CwSW ’ 
mlgw^w ~ s w) 

32n 2 rriy V cw 


9w9z 


1287r 2 m^,cvc ’ 
m T em z r]^ P3 
4 tt 2 °w s w 


cos(a — (3), 


(40) 

(41) 

(42) 

(43) 


where Cw, sw , tow, mz are the cosine, sine of the Weinberg angle, W and Z gauge boson mass, respectively. We 
omit A e because the Feynman diagram (e) gives no contribution to weak magnetic dipole moment. The Gf lP2P3 (x, y ) 
functions are: 


where 


Ga hh (x,y) = -^(C v {x + y)(mi i {x + y~l)-m h ))/M^, 

Gb ±H±u { x i y) = Am liX (x-l)/Ml 

Gd ±vv {x, V ) = —8 m h (x + y- l)(a; + y)/M%, 

G z f HT (x,y) =- [((x + y- l) 2 {x + y- 1) 

J m z z Mf 

+ Mf(log(Mf)(2- 3(x + y)) -x-y- 1) + m 2 z x(y{x + y - 1) +2)), 


M a = rnf.(x + y-l)(x + y)+mi j (x + y)-m 2 z xy-m 2 H (x + y-l) 7 
M b = m 2 .(x-l)x + m 2 z y{x + y-l)-m 2 H (x-l), 

M d = mf.(x + y-l)(x + y)-m 2 z xy-m 2 H (x + y-l), 

Mj = mf.((x + y-l)(x + y-l))-m 2 z x(y-l) + m 2 H y. 


The total contribution to (a ™) is: 


( fl ff)TflDM-//7 


(W\H ± H ± v 
\ a r ) 


£ (a?)* 1 * 

H=h°,H°,A° 


+ (d. 


W'jH^ 


{a. 


W^ZHt 


(44) 

(45) 

(46) 

(47) 


(48) 

(49) 

(50) 

(51) 


(52) 
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IV. CONSTRAINTS FROM LOW ENERGY PROCESSES 

Due to the introduction of a second scalar doublet, there are additional Feynman diagrams that are not present in the 
SM, involving the exchange of Higgs bosons ( h°, H°, A 0 , H ± ). There are currently strict bounds and measurements 
that have been obtained through various low energy phenomenological studies, this allows us to constraint parameters 
of the model such as tan fj(= tp) = ty j'U\ which controls the intensity of the coupling of the Higgs bosons to the 
fermions. In this section we show the allowed parameters space in the plane tp — the other parameters are 

fixed. In particular for a and (3 angles we study the scenario where a — /3 = ^, which reproduces the case where 
the coupling of h° to fermions is SM-like. In general, the favored region for THDM’s are phenomenologically viable 
for the scenario a — f3 = ^ m- We consider the relevant low energy constraints and the muon anomalous magnetic 
dipole moment for constraint the model parameters and then we give analytical predictions to magnetic and weak 
magnetic dipole moments of charged leptons in the sectionjv] using values of the allowed region. For the low energy 
constraints mediated by neutral Higgs bosons, we consider B® —> pp, p —> ey, r —> /xy, r —> ey, K — K mixing and 
h —■► Ijlkh and the rare Higgs decay h —>• Tp,. While for processes mediated by charged Higgs bosons b —» sy and 
B —» D(D*)tv. The processes H° —» 2 jets, Z —► //, e + e~ —> H°Z, H°A°, e + e~ -A bbH impose a limit on the mass 
of the Higgs boson given by m^o > 110.6 GeV at 95% [2E| and in the range 1 — 55 GeV is excluded [25] , 

On the other hand, a theoretical approach that can be used to try and understand the underlying physics is through 
the effective Hamiltonian, 


Weff — Vc KM Ci(p)Qi, (53) 

where Qi are the operators which govern the decay and Ci(p) are Wilson Coefficients. The scale p separates the 
physics contributions into short distance contributions contained in Ci(p) while the long distance contributions are 
contained in Qi but the Wilson Coefficients depend on the couplings of the particular model, which in our case are 
functions of a, /3, Qi, m,H(H ± ) an( i If - This will allow us to constrain the model parameters. 


A. Low energy constraints mediated by neutral Higgs bosons 

l. Bg —I PP 

A measurement of B® —> pp was reported in 135] . whose value is 

BR(B° S -a pp) = (2.8±°;*) x 10” 9 , (54) 

while the SM value is given by 136 

BR(B° -A pp) S M = (3.23 ± 0.27) x 10" 9 , (55) 

which puts very stringent constraints on the model. The Feynman diagram for this process at the quark level is shown 
in FIG. (3). 

S \ 

h xS 



Figure 3: Feynman diagram for the H 9 —► pp process at the quark level. The labels indicate the new contributions 
coming from the THDM-III. We will omit these labels in all further Feynman diagrams. 

We observe that the new contributions from Higgs bosons appear in the Hsb and Hpp couplings. The Branching 
ratio for this decay is given by E3, 
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BR ( B° s -A fi») = G f^M B r B \V tb V*\ 2 
2 ( 4m'/. \ 

' F «'{'-Ml +If>l 


\ x _ 4 ml 


M 2 B 


where F B and F^o are the form factors given by, 


(56) 


F h .ao = ~AM 2 B f Bs - ^C HA o, 

2 TYib + m s 

(57) 

where Gf is the Fermi constant, t b the lifetime of the B meson, M B the B meson mass, f Bs = 0.256 and C B , Ca° 
the Wilson coefficients that appear in the effective Hamiltonian, 


H eff = -2V2G F V tb V* Y, CiQu 

i 

where Qi are the effective operators and are given by: 

e 2 e 2 

Qh = j^mb{sPnb)(fln), Q a o = j^m b (sP R b)(p,j 5 n) 

and the Wilson coefficients C, are, 


(58) 


Ch = 


27TTO„ 


VtZVtbOe 


4 m 2 Tj 0 — 

H H=H°,h° 


E 


h 

lbs I mm 


Ca 0 = 


2tt m u 


V t * s V tb a em 


4m ^ 0 


Vbs Vmm 


(59) 


where m b , V * b ^ ts -,, m B o , myp , TOa° are bottom quark mass, CKM elements, heavy, SM-like, pseudo scalar Higgs bosons 

H,A° i 


mass, respectively. We can see through 


59 


the dependence on the model parameters since r] 


and therefore C b ,a° = C b ,a°(ol-, /3, $i, m#, 7 /). 


H,A° 

ff 


= V 


ff 


(a, /3, Oi, 7 /) 


2. Radiative decays U -A lj 7 

The second processes that we considered are the radiative decays li -A Ijj. The Feynman diagram with the largest 
contributions for this process is shown in FIG. [I], 


k 



h 


Figure 4: Feynman diagram for the h -A Ij'y process. 


The branching ratio for the general decay of a lepton (Z,) to a lepton of a different family (lj) is given by: j35[ 
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Br{k ->• Ij'y) = 


47rr i 


(\c 1 r '\ 2 + \c l [ li \ 2 ) , 


where the Wilson Coefficients C l ^' L are: 


°RH — 


LRH* LRH , LRH* LRH 
VliL VliL -r Vlili HljU 




such that H = h°,H°,A°. Cp is obtained by simply interchanging R by L. 


r 8 

r 7 

,-13 


The experimental bounds for each of the different processes are 

BR(t —> /zy) < 4.5 x 10 
BR{t —> ey) < 1.1 x 10 
BR(fi —y ey) < 5.7 x 10 

Again we can see that = C^(a, /3, 0j, mjj,y/). 

5. K — K mixing 

The amplitude for KK Mixing receives contributions from the Feynman diagram shown in FIG. [ 5 ]. 


(60) 


(61) 





Figure 5: Feynman diagram for the KK — mixing. 


The effective hamiltonian for this process at the quark level is written from [32] as follows 


H. 


AS=2 

eff 


G\M^ 

167T 2 




with (i = L , R , LR), where 


of 85 = (sP L d)(sP L d), 

Q sd,sd = ^p R( ^^p Rd ^ 

Q-sdM = ^p Ld j^p Rd y 


(62) 


(63) 

(64) 

(65) 
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and the Wilson coefficients are written here for a general multi-Higgs doublet model, 

c sd,sd _ 1671-2 \ H\2 y' (^la) 

/-<2 2 m s m d|»Ul 2 ’ 

'~ r F ,n W a=1 "‘if 


_ 1671-2 I |2 Y' ( Ula ) 

— r 2 m 2 msmd \^sb\ 2 

' J F rn W a=1 


C Lfl Sd = -7^Y msmd l^| 2 

^ F ,n W 


ULUla 


a—1 


°H 


( 66 ) 

(67) 

( 68 ) 


where i7 = i?°, /i°, H°. Ub a is the rotation matrix to go from weak to mass-eigenstates in the neutral CP-even Higgs 
sector, since we are working in general THDM, the rotation matrix for neutral Higgs bosons takes the form 


U = 


cos a — sin a 
sin a cos a 


Flavor violation arises from the non-diagonal terms of the K (neutral kaon) mass matrix, in particular the component 
My 2 ■ whose experimental value has been measured to be 


m k = ^Mk_ = 7 2948 x io- 15 GeF, 


but given that A Mk is obtained through 


M k 


(69) 


A M k = 2Re{K' } \H*/ f = 2 \K u ) 

C 2 A/t 2 

- f w -M k F 2 kV2 B k 


12tt 2 

t/- isd.sd , r) (/-isd.sd , /-isd.sd 
F 2 ,lrC lr + Hi t L [C L + C R 


(70) 


where Fk = 160 MeV, Mk = 497.6 MeV, ry 2 = 0.57, Bk = 0.85 ± 0.15, P 2 ,l_r = 30.6 and P\ j, = —9.3., we have 


My 2 = -FKT]iBK{m d m s )- 


E 

a—l 


(vcospy 

F> U 2 aUla jj 

H 2 ,Li?-o-b Hi L 

m H 


TT 2 TT 2 

Ur >n . 


+ 


L H 


l H 


(7!) 


Jf.. li y Ijlklk 

The Feynman diagram for the process /, —> Ijlklk is shown in the FIG.[6j, 


k h 

■*-^-*— 


H 


Figure 6: Feynman diagram for the li —> Ijlklk process. 
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The Branching ratio for this process including contributions from the three Higgs bosons is given by |43] 


Br{k 


Ijlklk') — 


5S h i k +2 Ti rn^rn^mf. (cos 2 (a-/3)sin 2 a 


2 11 tt 3 


TO 


h° 


sin 2 (a — /3)cos 2 a cos(a — (3)sin(a — cosa sina 


- 2 


l H° 




sin 2 \ \rujf_ 
m \o J 2cos 4 /? 

where r, is the time life of the U particle. The current upper bounds from PDG [44] are: 


BR(t —)• 
BR(/j, —» eee) 


< 1.0 x 10 -12 

< 2.1 x 10” 8 


B. Low energy constraints mediated by Charged Higgs bosons 

1 . b —» sy 

In the Standard Model the dominant decay of b quarks are given via couplings to the W~. And rare decays like 
b —>• sy are described by penguin diagrams involving W~ at one loop level, but in the THDM heavy charged Higgses 
also contribute at the one loop-level so the decay b —> sj remains as a good candidate to probe BSM physics. 

The Feynman diagram for b —> sy involving H~ is shown in FIG. [7]. 



Figure 7: Feynman diagram for the b —> S'y. 


The decay width is given by 


r( &-> si) = ^\V t :V tb \ 2 a em m 5 b \D\ 2 (72) 

where the reduced amplitude D and the Wilson coefficients are presented in [55], and discussed in Rt)j . 


2. B -> D(D*)tv 

The BABAR collaboration m measured the branching ratios of the B decays BR(B —> D(D*)tv ), finding: 

R(D) = 

R(D*) = 


0.44 ± 0.058 ± 0.042 
0.332 ± 0.024 ±0.018. 


( 73 ) 











where 
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R(D(D*)) 


BR(B -> D(D*)tv) 
BR(B ->• D{D*)lv) ’ 


with l = e, n, while the SM theoretical prediction for these processes is given by 


R sm (D) = 0.297 ±0.07, 

R sm {D *) = 0.252 ± 0.003 ±0.018, 


(74) 


(75) 


taken together the results give a discrepancy of 3.4 ct. These processes are present in other models that involve charged 
scalar bosons and can be used to constrain the parameters of such models, like in our case of the THDM-III. 

In this model the Feynman diagram that contributes to this process is shown in FIG. [8] 



Figure 8: Feynman diagram for the B —> D(D*)tv process. 


These ratios can be expressed as follows, 


R(D) = Rsm(D) 

x [ 1 ± 1.5 Re 


\Cr+C c l] 

-i-1 n 

ct + c? 

r^cb 

L U 5M J 


r^cb 
°SM 


R(D*) = R S m{D*) 

x ( 1 ± 0.12 Re 


cb 


C% - C' L 

r^cb 
°SM 


±0.05 


C cb sicb ^ \ 

R ~ 


/"icb 


(76) 


(77) 


where Cf R are the Wilson coefficients given in [48] . 

We find that this is the single most stringent low energy processes. Nonetheless, we are able to satisfy both R{D) 
and R(D*) within our model. The combined allowed region of charged Higgs boson mass (mj) and tp for all the 
previous low energy processes are shown in the following section. 


C. h° ->• r/i 

The flavour violating decay h° —>• r/i currently has an excess of signal events with a significance of 2.4 a [4j3], whose 
branching ratio reported by CMS with a 95% confidence level is given by: 

BR(h° ->• r/i) < 1.51%. 

The tree level Feynman diagram for this process is shown in FIG. [9]. 


( 78 ) 
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(a) 

Figure 9: Feynman diagram of hr — > Tfi. 


It is important to determine if our model acomodates this anomaly, to this purpuse we calculated the decay width 
given by: 


T(h° —> Tfj.) = 


N r 


8Tnn 2 h0 


{m 2 h o - (m r + nifj ,) 2 ) 3 ^ 2 
x ( m h° ~ (m T - m M ) 2 ) 1/2 , 


(79) 


where N c = 1 is the color number and m^o, m T , m^ are the neutral light Higgs boson mass, tau mass and muon mass, 
respectively. We found that the Branching ratio for this process as a function of tp, falls within the LHC reported 
measurement as shown in FIG. 10 , 



(a) 

Figure 10: Branching ratio of h° —» r/i decay as a function of tp. 

It might seem that the branching ratio for case 2a could grow in an undesirable way, but it becomes stable at the 
value of 2 x 10 —2 , while for the other two cases it decreases as tp increases. The allowed range over tp that we find in 
this process matches the allowed range obtained from the low energy process that will be presented in the following 
section. 


D. Allowed regions 


tp and mu 

We determine the allo wed parameter space by first finding the allowed region for all of the low energy process: |IV A 1[ 


IV A 2 
worked' 


IV A 3 


IV A 4 


and we then intersect this region with the allowed space of . In the analysis we have 

with the 7 ; = 0.5, {a — /3) = | scenario, mjj± = 500 GeV ,mn = TOjjc^m^o, 6 i^ = 0 and m^o = 125 GeV; 
the selection of these values are reasonable, since they satisfy current experimental bounds [33] on both neutral and 
charged heavy Higgs. 
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The allowed regions ( tpvsmn) for 


case 1, 2a, 2b is presented in FIG. 


11 



Figure 11 : Allowed regions for the cases: (a) case 1, (b) case 2a, (c) case 2b. 


Case a) has the largest allowed region, since it allows for small Higgs masses with a small tp ~ 10 through heavy 
Higgs of a large tp ss 50. While for case 2a we have a much more restricted parameter space, since small masses are 
allowed for a small range of tp ~ 6 — 8 and the heaviest higgses are between tp a, 12 — 16. For the case 2b we also 
see a small allowed zone of parameter space, except that unlike the previous figure, slightly larger values of tp are 
allowed. 


tp and m H ± 


For the charged heavy Higgs bosons we use the constraints |IV B ifTV B 2| to determine the allowed parameter space 
tpvsm H ±. Where we have used 7 = 0.5, (a — /3) = f = 500 GeV, m^ = 125 GeV and = §, #2 = f • The 


allowed regions ( tpvsm H ±) for case 1, 2a, 2b is presented in FIG. 12 



Constraints on the THDM must also satisfy electroweak precision data (EWPD). Thus, for our results to be reliable 
we verified that the parameter space regions explored are consistent with EWPD. In particular we found Higgs 
boson masses and tp that satisfy radiative corrections to the W boson mass in the range of 100 < rrifj < 900 GeV 
{mu = nine = m#) and 1 < tp < 40 for different values of m H ± = 300, 500, 800 GeV and mw = 125 GeV. Our 


results agree with the authors of f50l ■ A brief discution is included in App VIIB 
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V. RESULTS 

In this section we show the results obtained for the magnetic and weak magnetic dipole moments in the framework of 
the THDM-III with four zero textures. In the previous section we presented the low energy processes involving either 
charged or neutral Higgs bosons and we explored what are the available allowed regions that satisfy the experimental 
bounds or measurements. 


A. Muon anomalous magnetic dipole moment. 


In FIG. 13 we show the contributions to a M coming from the one loop diagrams in the context of THDM-III. 





(a) (b) (c) 

Figure 13: New contribution of the THDM-III to as function of Higgs mass (m#) and different values of tp for 
reasonable values of free parameters of the model: 7 = 0.5, (a — /3) = ^ , to # = to#o = 171^0, rrijj± = 500 GeV and 
To; jo = 125 GeV for: (a) case 1, (b) case 2a, (c) case 2b. The parallel lines in all the plots correspond to allowed interval 
of the discrepancy Abetween the prediction of the SM and experimental measurements with a 95% confidence 
level. 


In each of the plots of FIG. 13 we notice that small masses mu are allowed for small values of tp, while heavier 
Higgs masses are only allowed for greater values of tp. 


B. Magnetic and weak magnetic dipole moments of the r. 


Using the previously obtained parameter space, we show in FIG. 14 and FIG. 15 , the THDM-III scalar contribu¬ 
tions to the magnetic and weak magnetic dipole moments of the r. 


2x1 O' 7 
1x1 O' 7 
1x1 O' 7 
1x1 O' 7 
1x1 O' 7 
nf 1x1 O' 7 
9x1 O' 8 
8x1 O' 8 
7x1 O' 8 
6x1 O' 8 

sxio -8 ,,,,,,,,, 
200 300 400 500 600 700 800 900 1000 

m H (GeV) 



5x1 O' 8 
5x1 O' 8 
4x10‘ 8 
4x1 O' 8 
kT3x10' 8 
2x1 O' 8 
2x1 O' 8 
2x1 O' 8 

IxlO" 8 ■ ,,,,,,,,, 

200 300 400 500 600 700 800 900 1000 

m H (GeV) 



6x10 
6x10 
5x10 
5x10 
4x10 
of 4x10 
3x10 
2x10 
2x10 
2x10 
1x10 . 

200 300 400 500 600 700 800 900 1000 

m H (GeV) 



(a) 


(b) 


(c) 


Figure 14: Plots of the scalar contributions to a T in the THDM-III, (a) case 1, (b) case 2a, (c) case 2b. 
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from FIG. 14 , we find that the contribution to a T is contained between the current experimental value of [5], 

- 0.052 < a T < 0.058. (80) 

We can see that the order of magnitude varies moderately for each of our cases. The contributions of case 1 are within 
the range 0(1O -8 — 10~ 1 2 3 4 5 6 7 ) while case 2a and 2b are C?(10 -8 ), 0(1O -7 ), respectively. It is understandable because the 
three cases give similar contri but ions due to the structure of the Yukawa matrices we have considered. 

On the other hand, in FIG. 15 shown planes a^' vs run 


1x1 O' 8 - 

. 'r 20 x 

1x1 O' 8 - y 1^=30 * 

tp=40 , 

8x1 O' 9 - * • 

6x1 O' 9 - 

X • 

t __ x • 

® 4x10" 9 - X * • • 

2x1 O' 9 - * 

0x10° -1-1-1-1-1-1-1- 

200 300 400 500 600 700 800 900 1000 

m H (GeV) 


3x10‘ 9 
2x10‘ 9 
2x10‘ 9 
2x1 O' 9 
2x1 O' 9 
2x1 O' 9 
Sj, 1x1 O' 9 
1x10‘ 9 
IxlO" 9 
8x10' 10 
6x1 O' 10 

200 300 400 500 600 700 800 900 1000 

m H (GeV) 


2x10" 7 - * • 

2x10" 7 - * 

5 nf 2x10" 7 - ^ 

IxK)' 7 - X x x * ^ * 

5x1 O' 8 - * * 

0 x 10 ° -I-1-1-1-1-1-1-1- 

200 300 400 500 600 700 800 900 1000 

m H (GeV) 


(a) 


(b) 


(c) 



Figure 15: Plots of the scalar contributions to a^ in the THDM-III, (a) case 1, (b) case 2a, (c) case 2b. 

We note that the contributions of case 1 are 0(1O~ 9 — 10 -8 ) while case 2a and 2b are O(l0~ 10 — 10 -9 ) — 0(1O -7 ), 
respectively. 


VI. CONCLUSIONS 

We found that one loop level scalar contributions from the THDM-III with four zero textures can be a possible origin 
of the discrepancy A a^. For small values of the mass of the heavy Higgs run, tp is also small, and reciprocally, larger 
masses of run are permissible only for large tp. Furthermore, using current low energy constraints and a^ HDAI ^ ln , 
we obtained the allowed parameter space of the Higgs boson masses and tp which are also consistent with electroweak 
constraints. We determined that the magnetic dipole moment of the r lepton is of the order of 0(1O” 8 — 10” 7 ) and 
the weak magnetic of a^ is of the order 0(10” 10 — 10”'). Of the three cases that we consider the dipole moments 
found in case 2b have the largest value for a[ W \ 


Acknowledgments 

We acknowledge support from CONACYT (Mexico). And appreciate the advice and discussions of this work with 
Gilberto Tavares-Velasco, Lorenzo Diaz-Cruz and Agustin Moyotl. 


[1] G. Bennett et al. (Muon G-2 Collaboration), Phys.Rev.D73, 072003 (2006), arXiv:hep-ex/0602035. 

[2] M. Jung and A. Pich, JHEP1404, 076 (2014) [arXiv: 1308.6283 [hep-ph]. 

[3] D. Gomez-Dumm and G. Gonzalez-Sprinberg, Eur.Phys.J.Cll, 293 (1999), arXiv:hep-ph/9905213 [hep-ph], 

[4] J. Abdallah et al. [DELPHI Collaboration], Eur. Phys. J. C 35 (2004) 159. 

[5] S. Eidelman and M. Passera, Mod.Phys.Lett. A22, 159 (2007), arXiv:hep-ph/0701260 [hep-ph]. 

[6] Bolanos, A. and Moyotl, A. and Tavares-Velasco, G., Phys. Rev. D89 2014, arxiv: 1312.6860 

[7] T. Ibrahim and P. Nath, Phys.Rev. D78 , 075013 (2008), arXiv:0806.3880 [hep-ph], 

[8] A. Moyotl and G. Tavares-Velasco, Phys.Rev.D86, 013014 (2012), arXiv:1210.1994 [hep-ph] 

[9] J. Bernabeu, G. Gonzalez-Sprinberg, M. Tung, and J. Vidal, Nucl.Phys.B436, 474 (1995), arXiv:hep-ph/9411289 [hep-ph] 

[10] B. de Carlos and J. Moreno, Nucl.Phys. B519 , 101 (1998), arXiv:hep-ph/9707487 [hep-ph]; W. Hollik, J. I. Illana, S. 

Rigolin, and D. Stockinger, Phys.Lett. B416 , 345 (1998), arXiv:hep-ph/9707437 [hep-ph]. 










19 


[11] A. Heister et al. (ALEPH Collaboration), Eur.Phys.J. C30, 291 (2003), arXiv:hep-ex/0209066. 

[12] M. Laursen, M. A. Samuel, and A. Sen, Phys.Rev.D29, 2652 (1984). 

[13] Qing-Jun Xu, Chao-Hsi Chang, arXiv:1410.1774v2 [hep-ph] 

[14] Buras, Andrzej J. and Carlucci, Maria Valentina and Gori, Stefania and Isidori, Gino. JHEP, 10, 2010. 

[15] T. P. Cheng and M. Slier, Phys. Rev. D 35, 3484 (1987). 

[16] Honda, Mizue and Kaneko, Satoru and Tanimoto, Morimitsu, Phys. Lett. B593. 2004. 

[17] Branco, G. C. and Emmanuel-Costa, D. and Simoes, C. Phys. Lett., B690, 2010. 

[18] G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Rcbelo, M. Slier, and J. P. Silva, Theory and phenomenology of two- 
Higgs-doublet models, Phys. Rcpt. 516(2012) 1 [arXiv: 1106.0034 [hep-ph]]. 

[19] H. Nishiura, K. Matsuda and T. Fukuyama, PHYSICAL REVIEW D, VOLUME 60, 013006. 

[20] J.L. Diaz-Cruz, R. Noriega-Papaqui and A. Rosado, 

[21] G. C. Branco, D. Emmanuel-Costa and R. Gonzalez Felipe, Phys. Lett. B 477, 147 (2000) [hep-ph/9911418]. 

[22] M. Gomez-Bock and R. Noriega-Papaqui, J. Phys. G 32, 761 (2006) [hep-ph/0509353]. 

[23] J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, “The Higgs Hunter’s Guide,” Front. Phys. 80, 1 (2000). 

[24] D. Chang, W.-F. Chang, C.-H. Chou, and W.-Y. Keung, Phys.Rev. D 63, 091301(R)(2001). 

[25] F. Larios, G- Tavares-Velasco, C.-P Yuan, Phys-Rev. D 64, 055004. 

[26] OPAL Collaboration, G. Abbiendi et.al., Fur. Phys. J. C 18, 425 (2001). 

[27] Dean Carmi, Adam Falkowski, Eric Kuflik, Tomer Volansky, Jure Zupan, 10.1007/JHEP10(2012)196, arXiv: 1207.1718v3 
[hep-ph], 

[28] J. Abdallah et al. [DELPHI Collaboration], Eur. Phys. J. C 38 (2004) 1 [hep-ex/0410017]. 

[29] The OPAL Collaboration, Search for a low mass CP-odd Higgs boson in electron positron collisions with the OPAL detector 
at LEP2, Eur. Phys. J. C27 (2003), [hep-ex/0209068]. 

[30] M. Arroyo, J. L. Diaz-Cruz, E. Diaz and J. A Orduz-Ducuara, arXiv: 1306.2343 [hep-ph] 

[31] G. CVETIC, S. S. HWANG, and C. S. KIM, Int. J. Mod. Phys. A 14, 769 (1999). DOI: 10.1142/S0217751X99000385 

[32] Zhi-zhong Xing, Zhen-hua Zhao. D01:10.1016/j.nuclphysb.2015.05.027. arXiv:1501.06346v2 [hep-ph]. 

[33] P.M. Ferreira, R. Guedes, J.F. Gunion, H.E. Haber, M.O.P. Sampaio, R. Santos, The Wrong Sign limit in the 2HDM, 
arXiv:1410.1926 [hep-ph]. 

[34] A. Heister et al. (ALEPH Collaboration), Eur.Phys.J. C30, 291 (2003), arXiv:hep-ex/0209066 [hep-ex]. 

[35] V. Khachatryan et al., [CMS and LHCb Collaborations], arXiv:1411.4413 [hep-ex]. 

[36] C. Bobetli et al. , Phys. Rev. Lett. 112 (2014) 101801. 

[37] A. Dedes and A. Pilaftsis, Phys. Rev. D 67, 015012 (2003) [hep-ph/0209306]. 

[38] Crivellin A Kokulu A and Greub C 2013 Flavor-phenomenology of two-Higgs-doublet models with generic Yukawa struc¬ 
ture,” Phys. Rev.D879, 094031 [arXiv: 1303.5877 [hep-ph]] 

[39] J. Adam et al. [MEG Collaboration], Phys. Rev. Lett. 110, no. 20, 201801 (2013) [arXiv: 1303.0754 [hep-ex]]. 

[40] B. Aubert et al. [BaBar Collaboration], Phys. Rev. Lett. 104, 021802 (2010) [arXiv:0908.2381 [hep-ex]]. 

[41] K. Hayasaka et al. [Belle Collaboration], Phys. Lett. B 666, 16 (2008) [arXiv:0705.0650 [hep-ex]]. 

[42] Dedes, Athanasios and Pilaftsis, Apostolos, Phys. Rev., D67, 2003, 015012, 10.1103/PhysRevD.67.015012. 

[43] Y-F. Zhou, J.Phys.G30:783-792,2004. 

[44] Olive, K. A. and others, Review of Particle Physics, Particle Data Group, Chin. Phys., C38, 2014. 

[45] F. Borzumati and C. Greub, Phys. Rev. D 58, 074004 (1998) [liep-ph/9802391]. 

[46] P. Ciafaloni, A. Romanino and A. Strumia, Nucl. Phys. B 524, 361 (1998) [hep-ph/9710312]. 

[47] Lees, J. P. et al. Phys. Rev.D88, 2013, 7, 072012. 

[48] A. Crivellin, C. Greub and A. Kokulu, Phys. Rev. D 86, 054014 (2012) [arXiv: 1206.2634 [hep-ph]. 

[49] V. Khachatryan et al. [CMS Collaboration], arXiv:1502.07400 [hep-ex]. 

[50] Alessandro Broggio, Eung Jin Chun, Massimo Passera, Ketan M. Patel, Sudhir K. Vempati, arXiv: 1409.3199vl [hep-ph]. 

[51] David Lopez-Val and Joan Sola, arXiv:1211.0311v3. 



VII. APPENDIX 


A. Yukawa elements for Case 1 


>21 = 


0 — i(0l+02) 


c*(01+02) 


1 7mlra 2 (77712 + m 2 — 7773) 

(mi + m 2 ) (mi + m 3 ) (7 to 2 — to 3 ) 
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(mi + m 2 ) (m 2 - m 3 ) (77772 - m 3 ) 
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^23 = 
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= -i(e 1 +e 2 ) / m 2 ( 7 m 2 + m 2 - m 3 ) / mim 2 mj 


(mi + m 2 ) (m 2 - m 3 ) y 7 m 2 - m 3 


7 mim 2 


(m 2 — m 3 ) (mi + m 3 ) ( 7 m 2 — m 3 ) 
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Y 3 1 = e ~ i{ei+02) 


x i - 


7TO1TO2 (7m 2 + ?n 2 - m 3 ) 

(mi + m 2 ) (mi + m 3 ) (7m 2 — m 3 ) 

m 3 (7m 2 + ?n 2 - m 3 ) (—7m 2 + mi + m 3 ) 


i( 0 i+ 0 2 ) 


+ e' 




(m 2 — m 3 ) (mi + m 3 ) (7m 2 — m 3 ) 
7m 2 ?n 3 


(m 3 - 7m 2 ) 


(m 2 - m 3 ) (mi + m 3 ) 


/7m 2 (7m 2 + m 2 — m 3 ) (—7m 2 + mi + m 3 ) 
7m 2 — m 3 


/mg (~7m 2 +mi + m 3 ) / i(9i0: 


—1— g 

V (mi + m 2 ) (mi + m 3 ) 
x ( 7 m 2 — mi + m 2 ) — e *^ 201+ei+02 ^ 

/ 7mim^ 

V (m 2 - m 3 ) (mi + m 3 ) (7m 2 - m 3 ) 


7 m 2 m 3 


(m 2 — m 3 ) (mi + m 3 ) 


mim 2 m 3 
7m 2 — m 3 


„i(ei+e 2 ) h m 2 (7^2 + m2 - m 3 ) (~7m 2 + mi + m 3 ) 
6 1 / 


7 m 2 — m 3 


x 1 - 


m 3 (7m 2 + ?n 2 — m 3 ) (—7m. 2 + mi + m 3 ) 


+ e 


(m 2 — m 3 ) (mi + m 3 ) (7m 2 — m 3 ) 
_ i( e 1+ e 2 ) / mim 2 m 3 / 7 m 2 m 3 


7m 2 - m 3 Y (m 2 - m 3 ) (mi + m 3 ) 
m 2 m 3 (7m 2 — mi — m 3 ) 


(mi + m 2 ) (mi + m 3 ) ( 7 m 2 - m 3 ) 
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+32 = 


X W- 


+ e 


-e i6l62 \ - 

V ( 

7 m| (— 7 m 2 + mi + m 3 ) ( 

[mi + m 2 ) (m 2 - m 3 ) (71712 - m 3 ) V 

/ m 3 (7777-2 + m 2 - m 3 ) (—7777-2 + mi + m 3 ) 

/ (?T 

— 2 ( 01 + 02 ) 

7-2 - m 3 ) (mi + m 3 ) (7m 2 - m 3 ) 

/ 7 m 2 m 3 

V (777-2 - m 3 ) (mi + m 3 ) 

/7m 2 (7m 2 + m 2 - m 3 ) (—7m 2 + mi + m 3 ) 


iO-iQz 


(m 3 - 7 m 2 ) 


7 m 2 - m 3 


+ e -ie!fl 2 / m 2 ( 7 W 2 + m 2 - w 3 ) / ieie2 / 7 m 2 m 3 

V (mi + m-2) (m 2 - m 3 ) \ y ( m 2 “ m 3) K + m 3 ) 


x (7TO2 — mi + m 2 ) + e d++++++) — 


mim2m 3 
7m 2 - m 3 


7 ml mg 


(m2 - m 3 ) (mi + m 3 ) (7777-2 - m 3 ) 


+ e' 


;(fl 1 + e 2 ) / 7 to 2 (7rn 2 + m - 2 - m 3 ) (-77712 + m t + ?n 3 ) 

V 7m 2 - m 3 


X 4 - 


m 3 ( 7 m 2 + m 2 — m 3 ) (—7777-2 + mi + m 3 ) 


(m 2 — m 3 ) (mi + m 3 ) ( 7 m 2 - m 3 ) 


+ e ~ i{ei+e2 \- 


mim2m 3 


? 77 i m 3 (7m 2 + m 2 — m 3 ) 


x . / — 


7m 2 — m 3 y (mi + m 2 ) (m 2 - m 3 ) (7m 2 - m 3 ) 

7m 2 ?77 3 

(m 2 - m 3 ) (mi + m 3 ) 


+33 = 2 J- 


7 m 2 m 3 


(m 2 - m 3 ) (mi + m 3 ) 


mim2m 3 


77771)712 


7m 2 - m 3 

cos ( 0 i + 0 2 ) 


+ 


(m 2 - ? 77 - 3 ) (mi + ? 77 - 3 ) (7m 2 - 777 - 3 ) 

1 7m 2 (7m 2 + m 2 — m 3 ) (—77712 + mi + m 3 ) 
77712 - m 3 

m 3 (7m 2 + m-2 - m 3 ) (—7m 2 + mi + m 3 ) 


/ |«3 »2Tra2-»‘3 ">2traiTm3 , n . n n 1 

X 4 /- 7 - 77 - , w - 7 -- COS (-0201 + 01 + 02) 

y (m 2 - 777- 3 ) (mi + 777- 3 ) ( 7 m 2 - 777 - 3 ) J 

m 3 (m 3 (277712 — 777-1 + m 2 ) + m 2 (27ml — 27(7 + l)m 2 + mi) - m 3 ) 

(m 2 - ?? 7 3 ) (mi + m 3 ) 


+ 


B. Electroweak constraints 

Following the work done in [ 51 ] the W mass in the THDM is perturbatively calculated through the relation 

Ml ' 


M^ = 


1+4 1 - 


47ra<, 


V2 G f M 2 7 


(1 + Ar) , 


( 81 ) 
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where Gf , cc em , Mz are the Fermi constant, fine-structure constant and Z boson mass, respectively. The terms A r 
are radiative corrections to the W boson mass given by: 

A r = A r SM + A r THDM (82) 


where 

A r THDM = A a THDM - ^fSp THDM + A r rem , (83) 

s w 

such that Aa THDM is the photon vacuum polarization of the THDM, A r rem incorporates the remaining contributions 
and 


6p THDM = ( cos2 (ft - a) [F(ml o, m 2 H± ) - F(m 2 h0 , m \ 0 )] 

+ sin 2 (/3 - a) [F(m 2 H0 , m 2 H± ) - F(m 2 H0 , m 2 A0 )] + F(m 2 A0 , m 2 H ± ) 

- 3cos 2 (/3- a)[F(m,go,Mw) + F(m 2 h0 ,M 2 z )\ 

- F(m 2 H0 ,M 2 Z )~ F{m 2 h0 , Af^)) , (84) 


- if x ^ y (85) 

2 x - y 2/ 

F{x, y) = 0, if x = y. (86) 



